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Abstract 

We present a simple and non-invasive experimental procedure to measure the linear viscoelastic properties 
of cells by passive video particle tracking microrheology. In order to do this, a generalised Langevin 
equation is adopted to relate the time-dependent thermal fluctuations of a bead, chemically bound to the 
cell's exterior, to the frequency-dependent viscoelastic moduli of the cell. It is shown that these moduli 
are related to the cell's cytoskeletal structure, which in this work is changed by varying the solution 
osmolarity from iso- to hypo-osmotic conditions. At high frequencies, the viscoelastic moduli frequency 
dependence changes from cx uj^^^ found in iso-osmotic solutions to oc uj^/^ in hypo-osmotic solutions; the 
first situation is typical of bending modes in isotropic in vitro reconstituted F-actin networks, and the 
second could indicate that the restructured cytoskeleton behaves as a gel with dangling branches^\ The 
insights gained from this form of rheological analysis could prove to be a valuable addition to studies 
that address cellular physiology and pathology. 

Introduction 

As many authors have noted, the mechanical properties of a cell's cytoskeleton can influence factors 
such as growth, apoptosis, motility signal transduction and gene expression [1 . Related to this, there 
is a desire to be able to provide a rheological interpretation of the cell's viscoelastic response that has 
the potential to yield quantitative information on the cell's cytoskeletal structure. Consequently, in this 
work, we have developed a means of using passive video particle tracking microrheology measurements 
to quantitatively measure changes in the viscoelastic properties of a cell as a consequence (in this case) 
of simple changes in its external environment, i.e. subjecting a cell to a hypo-osmotic shock. 
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Osmotic regulation and the transport of osmotically active molecules are fundamental to both metabolic 
processes and homeostasis of cells and require precise regulation and maintenance of intracellular water 
[2]. The ability of many cells to regulate their volume, via internal restructuring, in response to osmotic 
changes in their environment is an essential component of normal cellular function. This regulatory vol- 
ume change is linked to a reorganisation of the cytoskeletal actin networks 01 1^ . Exposure to a hypotonic 
solution typically induces an initial rapid swelling followed by a shrinkage of the cell that occurs over 
a slower timescale of several minutes. The modulation of the actin dynamics and polymerisation that 
accompanies these regulatory volume changes has been shown in various cell types by methods such as 
DNase I inhibition assay and fluorescence measurements of phalloidin-labelled actin [3] . 

In this work we have developed a video particle-tracking tool to study the microrheology of cells, 
using a Jurkat lymphocyte cell line as a model system. In general, lymphocytes have been shown to have 
a response to hypo-osmotic solutions that is influenced by tyrosine kinase activity [5J and cytoskeleton 
(F-actin) participation in ion channel activation [6 . Furthermore, the interest in T-cells stems from 
their important role in the regulation of immune responses. In particular, Jurkat cells, as used in this 
study, are CD4-I- T-lymphoma cells that are often utilised as models to understand T-cell signalling [7] 
and HIV-l dissemination in viral pathogenesis |B]. 

The linear viscoelastic properties of a material can be represented by the frequency-dependent dy- 
namic complex modulus G* (lu) , which provides information on both the viscous and the elastic nature 
of the material (i.e. on how the matter dissipates and stores the energy transferred to it) at different 
frequencies uj; it is defined as the ratio between the Fourier transforms (denoted by the symbol """) of 
the applied stress a(t) (which is proportional to the applied force) and the resulting strain ^{t) (which is 
proportional to the material/cell deformation): 

^*^"^^^=/r7(t)e--dt 

where i is the imaginary unit (i.e. = ^1)- 

The standard method of measuring G*{uj) is based on the imposition of an oscillatory stress a{uj,t) 
given by a formula such as: 

a{uj,t) = (To sin(wt) (2) 
where ctq is the amplitude of the stress function, lo is the frequency and t the time. The resulting 
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oscillatory strain 7(w, i) then has the formula: 

7(cj, t) = 7o sin(wt + (5(a;)) (3) 

where 70 is the strain amplitude and S{uj) is the frequency-dependent phase shift between the stress and 
the strain. The amplitudes of the complex modulus in-phase and out-of-phase components are propor- 
tional to the ratio between amplitudes of the stress and of the strain, with constants of proportionality 
defining the storage (elastic) G'{uj) and the loss (viscous) G"{uj) moduh, respectively: 

G* (w) = G'(w) + iG"{uj) = ^ cos{5{uj)) + sin((5(w)) (4) 

70 70 

For example, in the case of a purely elastic solid, the stress and the strain arc in phase (i.e. Hooke's law, 
the material deformation is directly proportional to the applied force) and 5{uj) — 0; whereas, for a purely 
viscous fluid, such as water or glycerol, S{uj) ~ tt/2. For complex solids (e.g. gels, cells) or viscoelastic 
fluids (e.g. blood, saliva) S{uj) would take any value between these limits (i.e. ^ 5{uj) ^ 7r/2) depending 
on the frequency at which the force (stress) is applied. 

The aim of this article is to present a straightforward and non-invasive procedure for measuring the 
in vivo linear viscoelastic properties of single cells via passive video particle tracking microrheology of 
single beads attached to the cells' exterior. This method has advantages over both complicated active 
microrheology techniques, where complex experimental set-ups are necessary to exert an external force 
for performing stress-controlled measurements; and invasive passive video particle tracking of submicron- 
probes embedded (either via endocytosis or micropipette injection) within the cell's cytoskeleton (e.g. 
[HI [ini [TTJ [121 [131 [m inj • in particular, the procedure consists of measuring the thermal fluctuations of a 
bead chemically bound to the cell's exterior (Figure[T]), for a sufficiently long time. A generalised Langevin 
equation was adopted to relate the time-dependent bead trajectory, r(t), to the frequency-dependent 
moduli of the cell. Notably, the procedure presented here represents an alternative methodology that 
can be extended to many experimental formats and provides a simple addition to existing cellular phys- 
iology studies (e.g. those monitoring cell pharmacological response). Indeed, when compared to single 
cell viscoelasticity assays such as magnetic tweezers, atomic force microscopy and optical stretcher, as 
employed in [21 UHllIIl [HI [SI Ull [H] , our method has the advantage of revealing the changes of the cell's 
viscoelastic properties over a wide range of frequencies (here from ~0.6 Hz up to ~600 Hz), to a high 



4 



level of accuracy, whilst it experiences an induced physiological process. 




Figure 1. (Left) Top-view of a 5 iim diameter silica bead bound to the surface of a Jurkat 

cell. (Centre) Schematic side-view of the left image. (Right) The thermal fluctuations of a 5 iim silica 
bead chemically bound to the surface of a Jurkat cell. Images were analysed in real-time and the 
coordinates of the bead's centre of mass were stored directly on the random-access memory (RAM) of 
the computer, at time intervals down to l/600s, using our own suite of image analysis software written 
in Lab VIEW j23j. 



Analytical model 

The equation describing the pseudo Brownian fluctuation of the randomly varying bead position r(i)Vi 
can be derived by means of a generalised Langevin equation, which was originally derived to describe 
the Brownian motion of a freely diffusing particle suspended in a purely viscous fluid [16 and represents 
Newton's second law as written for stochastic processes: 

md{t) = fR{t)^ f [Ut - t) + Ut ~ T)]viT)dT (5) 
Jo 

where m is the mass of the particle, a{t) is its acceleration, v{t) is the bead velocity and /i?,(i) is the usual 
Gaussian white noise term that models stochastic thermal forces acting on the particle. The integral 
term represents the total damping force acting on the bead, which, based on the superposition principle, 
incorporates two generalised time-dependent memory functions Cc(i) and Cs(^) that are representative 
of the viscoelastic nature of the cell and the solvent, respectively. These memory functions are directly 
related to the materials' complex modulus as shown hereafter. 

In the case of the solvent, which in this work is the media surrounding the cell, the relationship 
between memory function and complex modulus is straightforward. Indeed, to a first approximation, 
using the assumptions adopted by Mason and Weitz [17] when studying the motion of thermally excited 
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free particles, at thermal equilibrium the complex viscoelastic modulus of the solvent is related to the 
memory function Cs(0 by the expression: G*{uj) = iujC,s{'^) /^t^R, where R is the bead radius and Qs{^) 
is the Fourier transform of Cs(^)- In the case of the cell, we assume a similar relationship between G* {uj) 
and Cc{t) to that given above for the solvent, but with a different constant of proportionality that we will 
call /3; i.e., G*(w) — iLoQci'-^)/ Note that, (3 may vary for different cells because it depends on (z) the 
cell radius (i?c), the number of the chemical bonds between the bead and the cell, {Hi) the contact 
area between the cell and the glass coverslip, and {iv) the relative position of the bead with respect to 
both the cell's equatorial plane and the glass coverslip. In this study we will only focus on the changes 
in cell dynamics due to variations in osmolarity. This can be achieved since all the above experimental 
parameters, with the exception of Rc, are (in good approximation) time-independent, and will not change 
significantly during the course of a set of measurements, as shown in Figure[4]for measurements performed 
on cells in iso-osmotic condition. In addition, it is important to highlight that the cells maintained a 
spherical shape well beyond the duration of the experiment. Moreover, the cells adhered to the coverslip 
and no drift was observed during the experiment. Finally, given that we measure an increase of the cell 
radius [Rc) of ^ 5% for the hypo-osmolarity experiments described below, we also make the assumption 
that this small change does not affect the dynamics of the system appreciably. 

We now describe how the thermal fluctuations of a bead, chemically bound to a cell, can be investigated 
to determine the viscoelastic properties of the cell through the analysis of the time dependence of the 
bead's mean-square displacement (MSD) <^Ar^(T)) = i\f{t + r) — r{t)f^ , in the radial direction of the 
cell; where t is the absolute time and r is the lag-time (i.e. time interval). The average is taken over 
all initial times t. Under these circumstances, it is an easy step to show that, at thermal equilibrium, 
Equation [5] can be reorganised to express the viscoelastic moduli of the cell as function of the Fourier 
transform of the mean-square displacement: 

' = — 1 (6) 

G'q iu)Il{uj) /3Gq /3Gg 

where Gq is the limiting value, for vanishingly low frequencies, of the elastic modulus of the whole system 
(i.e. cell plus solvent), which in this work is Gq = G*(a;) for w — > 0; and II(i^) is the Fourier transform 
of the normalised mean-square displacement n(r) = (Ar^(r)) / (r^), introduced by Tassieri et al. [TS] 
in the study of the thermal fluctuations of an optically trapped bead suspended in a viscoelastic fluid. 
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The term (r^) is the time-independent variance of r{t) measured for a sufficiently fong time. Moreover, 
it has been shown |18j that for a constrained bead (e.g. one that is optically trapped or, as in this work, 
chemically bonded to the cell) 11{t) = 1 at large time intervals. It is important to highlight that (as in 
the case of an optically trapped bead) the term /3Gq (representing the coefficient of the elastic restoring 
force) can be determined by means of the principle of equipartition of energy, which in one dimension is 
written as: 

ksT ^ f3G'o{r') (7) 

Usefully, two further simplifications can be made to Equation |6] because (a) for micron-sized silica 
beads, the inertia term mu!^ is negligible up to frequencies on the order of MHz and (6) for solvents having 
frequency-independent viscosity rj (e.g. water), G*(a;) simplifies to iujrjg. Thus, Equation |6] provides the 
viscoelastic properties of the cell (scaled by Gg) over a range of frequencies that is limited at the top end 
by the rate of acquisition of images to determine the bead position, and at the bottom end by a cutoff 
frequency given by LUcutoff — PG[^/riQ, where rjQ is the limiting value, for vanishing frequencies, of the 
viscosity of the whole system. 

Finally, in order to evaluate the Fourier transform in Equation [6] we adopt the analytical method 
introduced by Evans et al. [19] (i.e. Eq. (10) of Ref. [19]), which is applied directly to the experimental 
data points and has the advantage of removing the need for Laplace and inverse-Laplace transformations 
of experimental data [20] . 

Experimental details 

Anti-CD4+ (Invitrogen) was attached to 5/im carboxylate functionalised silica beads (Bangs Laboratories 
Inc) using a method based on a previously described protocol [21 . The adaptation of this method involved 
first binding a short polyethylene glycol chain to the beads so as to provide a flexible linker on to which 
the anti-CD4-|- could be bound. 

In detail, the microspheres were washed (3x) with 0.1 M Na2CO^ buffer (pH 9.6) using centrifugation 
between washes to remove the supernatant liquid. This ensured that all the carboxylate groups were 
deprotonated. Following the final wash, the beads were resuspended in RO water and centrifuged again. 
After removal of the water, the remaining pellet was resuspended in freshly prepared 0.02 M sodium phos- 
phate buffer (pH 6) containing 2 mM l-ethyl-3-(3-dimethylaminopropyl)-carbodiimide (EDC, Sigma) and 
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5 mM sulfo-N-hydroxysuccinimide (sulfo-NHS, Thermo-Fisher) . The microsphere solution was then ag- 
itated for 15 minutes to allow the sulfo-NHS to activate the carboxylate groups on the surface of the 
microspheres. The solution was centrifuged, then the supernatant removed and replaced by a 10 mM solu- 
tion of the flexible linker, 0-(2-Aminoethyl)-0'-(2-carboxyethyl)polyethylene glycol 3,000 hydrochloride 
(Aldrich) in pH 7.4 phosphate buffer. The activated microspheres were then agitated in the flexible linker 
solution for 2 h at 4°C, before further centrifugation and washing with pH 7.4 buffer. The anti-CD4-|- 
was then bound to the carboxylate terminals of the PEG linker from a 1 mg/ml solution of anti-CD4+, 
using the same EDC/sulfo-NHS activation procedure outlined above. 

The Jurkat cells were obtained from ATCC (clone E6-1, TIB-152) and maintained in RPMI 1640 
(Ix) media containing L-Glutamine and 25mM HEPES [Gibco, Invitrogen]. Prior to the experiment they 
were suspended in PBS [Gibco, Invitrogen] with an osmolarity of 275-304 mOsm, then mixed with the 
functionalised beads and allowed to equilibrate for ~10min. 

Beads were optically trapped and moved to the equatorial plane of the cell, where they were held for 
a sufRcient time (~2 min) until the binding process was completed. Trapping was achieved by using a 
CW Ti:sapphire laser system (M Squared, SolsTiS) which provides up to 1 W at 830 nm. Holographic 
optical traps were created via the use of a spatial light modulator (Boulder XY series) in the Fourier 
plane of the optical traps [52]. The microscope was equipped with a lOOx 1.3NA (Zeiss, Plan-Neofluor) 
objective lens, which was used both to focus the trapping laser beam and to visualise the particles via a 
Prosilica PCI-6713 camera. Samples were mounted on a motorised microscope stage (Prior Pro-Scan II). 
The beads' trajectories were measured and stored on a personal computer in real-time, at frequencies up 
to ~kHz (depending on the size of the analysed region of interest (ROI) within the camera field of view) , 
using our own suite of image analysis software written in Lab VIEW [53] . 

Results and Discussion 

In order to validate the analytical model introduced in this work, we measured the dynamic response of 
an analytical system that could be considered as an ideal viscoelastic medium with known storage and 
loss moduli; i.e. an optically trapped bead suspended in a Newtonian fluid. Indeed, for such a system, 
the elastic modulus must be proportional to the optical tweezers' (OT) trap stiffness (k), this being the 
only elastic component present in the system (i.e. the ensemble of Newtonian fluid, bead and OT); and 
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the viscous modulus must be linearly proportional to the frequency (w) with constant of proportionality 
given by the viscosity of the Newtonian fluid (77). 




Figure 2. (A) The mean-square displacements vs. lag-time of optically trapped beads in 
Newtonian fluids. The square and the triangle symbols refer to measurements performed with 2 /xm 
diameter silica bead in Glycerol/water mixtures at 10% and 20% w/w, having viscosity values of 1.15 
mPas and 1.54 mPas, and trap constants of 0.85 fiN/m and 0.81 ^N/m, respectively; The circle 
symbols refer to a measurement performed with 5 /xm diameter silica bead in water having viscosity 
value of 0.89 mPas and trap constants of 1.68 fiN/m (note that for this measurement only we used an 
oflt-line faster camera with internal memory; not suitable for continuous acquisition). The inset shows 
the same data as above but with the abscissa scaled by K/{6Trr]R). (B) The normalised viscoelastic 
moduli vs. the normalised frequency for all three analytical systems described above. 

In Figure [2J3 we present the normalised viscoelastic moduli versus the normalised frequency for three 
systems made up with diverse combinations of Newtonian fluids and beads having different viscosities 
and size, respectively. For each system, the moduli have been derived by the analysis of the thermal 
fluctuations of an optically trapped bead by means of Equation[6] wa the relative MSD shown in Figure[2]A_. 
In this case, the term G*{uj) in Equation |6] represents the viscoelastic nature of the analytical system; 
the term Gl{uj) simplifies to iajT] and is coincident with 3[G*(w)] (and therefore it does not need to be 
accounted for in Eq. |6| ; the term (3Gq = k is determined by means of the principle of equipartition of 
energy with /3 = 6ttR. In summary, for both moduli, the normalising factor is given by k/IGitR); whereas, 
for the abscissa, the scaling factor is given by a cut-off frequency deflned as k/ {QT:Rri), which also works 
in the time-domain as shown in the inset of Figure [2]^. In Figure [2|3 it is interesting to highlight that 
at low frequencies (i.e. up to the cut-off frequency) the elastic modulus G'{uj) is actually normalised 
to 1, as expected; whereas, at high frequencies (i.e. for lu > k/{QttR)) G'{uj) vanishes towards zero. The 
latter phenomenon is itself an interesting feature that is currently an object of discussion in literature 
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[24l [25l |26] . In particular, it is argued that the analytical models adopted so far for data analysis of 
the thermal fluctuation of an optically trapped bead suspended into a Newtonian fluid are unable to 
provide the high frequency elastic component of the optical tweezers (OT), because they would "suffer 
from systematic truncation errors at high frequencies". In contrast, we believe that the absence of the 
optical tweezers' high frequency elastic component is an intrinsic feature of the system made up by the 
ensemble of (Newtonian) fluid, bead and OT 27 ; on which we will provide a detailed explanation, but 
in a different work. On the other hand, the viscous modulus G"{lj) grows linearly with the frequency 
(thus in a Log-Log plot it is represented by a line with unity slope) over the entire experimental frequency 
window, with constant of proportionality given by the viscosity of the Newtonian fluid, as expected. Note 
that the results shown in Figure [2] are in good agreement with those published in literature [HJ [551 Ull 
for an optically trapped bead in water, but they used a different analytical method for data analysis and 
they overlooked the scaling process introduced in this work. 

Having demonstrated the validity of our analytical model, we now investigate its applicability to 
cellular physiology studies by focussing on changes in the cells' microrheology due to changes in osmolarity 
of the surrounding solution. 

It is anticipated that there are two major physiological processes that accompany the changes in 
osmolarity of the solution surrounding the Jurkat cell: (I) an actin cytoskeleton reorganisation and (II) 
an increase of the cell volume as the solution becomes hypo-osmotic. In order to quantify the changes 
in cell volume, we assumed that the cells were spherical, Figure [T] In particular. Figure [3] shows that 
the relative radius changes are < 6% when the surrounding PBS solution was made hypo-osmotic by 
addition of 10% v/v distilled water. The data of Figure [S] also shows that, after an initial swelling, the 
cells contract to a still swollen state that is ~5% larger than when they were in isotonic PBS. Based on 
these observations, which are in good agreement with those already existing in literature (e.g. |28j), we 
assume that these small radius changes do not affect the dynamics of the system appreciably. 

To study how the in vivo viscoelastic properties of Jurkat cells vary with the osmolarity of the solution, 
we have adopted the model described in the Analytical Section. In particular, analysis of different time 
regimes in the time dependent normalised mean-squared displacement, n(r), has the potential to reveal 
information both on the pure elastic component of the cell (at long time intervals or low frequencies) and 
on the fast dynamics occurring at small length scales (e.g. those related to the transverse bending modes 
of single actin filaments in the cytoskeleton). 
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Figure 3. Averaged relative cell radius changes vs. time as a consequence of the change in 
osmolarity of the cell solution; from iso-osmotic conditions in PBS buffer towards a 
hypo-osmotic condition obtained by adding 10% v/v of distilled water to the initial 
solution. The osmotic shock causes a rapid swelling of the cells, followed by a partial recovery of the 
cells' volume. These results are in good agreement with those reported in literature, e.g. |28j . 

Elastic properties of the cell derived from n(r) 

The first result that can be obtained through the analysis of the thermal fluctuations of the bound bead 
(Figure [T]) is the relative change of the cell's low frequency elastic plateau modulus G^(0), Figure [4] 
This can be evaluated from the time-independent variance (r^) of the constrained bead in the radial 
direction of the cell, by adopting the principle of equipartition of energy as in Equation [Tj In particular, 
although the absolute value of G^(0) can not be determined because of the unknown factor /3, the relative 
change between different values of /3G"^(0), obtained from n sequential measurements performed on the 
same cell, differ from each other only because of both the small increase of the cell radius (measured as 
~5%) and the variation of G^(0). Thus, from Figure |4] it is clear that a change of 10% in osmolarity 
towards hypotonicity, which occurs at t ~9min as a result of adding distilled water to the iso-osmotic 
solution, induces a temporary increase of /3G^(0) that lasts for approximately 15-30min. At its greatest, 
this change in the low frequency elastic modulus dependent term, /3G^(0), is over 3 fold in magnitude. 

The explanation for such a large increase may in part be due to the change in radius; however, if this 
were solely responsible, it would be expected that /3G^(0) would remain significantly above its original 
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Figure 4. The absolute value of /3G^(0) measured in isotonic (PBS) condition (circle 
symbols) and when the solution is made 10% hypotonic (square symbols) The arrow indicates 
the time at which the isotonic solution is made hypotonic. The dotted line is a guide for the eye, so 
that the later points in the graph can be compared to the first points measured in isotonic condition. 

value, as the cells do not regain their original radius. Thus, since /3G^(0) does return to close to the 
original value (at t ~25 min), consideration should be given to the temporary increase in electrostatic 
intermolecular interactions within the cell due to the reduction in ionic strength of the solution (i.e. 
a reduction in screening of the charges on cytoskeleton molecules by solution based ions). This latter 
phenomenon would also alter the cell's viscoelastic properties and any resulting actin cytoskeleton re- 
arrangement may be seen as part of an attempt of the cell to re-equilibrate the solution osmolarity, as 
discussed below. 

Evidence for an increase in intermolecular interactions in the cytoskeleton being the cause of the 
observed increase in /3G^(0) comes from studies of polymer gels [29l |30l [31] , for which the elastic plateau 
modulus (G'(0)) of an ideal cross-linked polymer network can be written: 



where kg is the Boltzmann constant, T is the absolute temperature, and -0 and / are the number 
of network strands and cross-links per unit volume, respectively (see Figure pi). Here, an increase in 



G'(0) ^ ipksT 



(/-2) 
/ 



(8) 
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number of strands and/or number of eross^inking interactions causes an overall increase in G'(0). Note 
that although it is appreciated that there may be more order associated with macromolecules within the 
cell's cytoskeleton, the generalised polymer network of Figure [5] nevertheless has many similarities with 
the complex structure of the cytoskeleton. For example, it is well known that all cells possess a host of 
actin-binding proteins that can align actin filaments into bundles and then cross-link filaments or bundles 
into networks |32j . Bundled actin and scruin in vitro networks act similarly to cross-linked networks and 
show a strong increase in elasticity as the degree of bundling increases [33j . 




Figure 5. Schematic view of a randomly cross-linked network with '•^dangling branches" 
defects. Circles represent cross-linking junctions and arrows denote attachments to the 
macroscopic network. The "dangling branches" do not contribute to the gel elasticity because they 
cannot bear stresses; in contrast, they can dissipate energy by friction with the solvent. 

Thus, by analogy with the models for a polymer network, a >3 fold increase in the cell elastic plateau 
modulus could be explained by an increase of intermolecular interactions (i.e. electrostatic interactions) 
that would have the effect of forming temporary cross-linkages and reshaping the isotropic cytoskeleton 
structure. This scenario is in good agreement with confocal microscopy results already reported for 
multiple cell lines and the observations that actin forms bundles in aortic bovine endothelial cells after 
hypotonic exposure |34) . 
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Fast dynamics of the cell derived from 11 (r) 

In order to study the high frequency behaviour of cells, we have used Equation |6j which directly relates 
the normalised mean-square displacement, TI{t), of the bead to the cell complex modulus, G*{uj). In 
particular, we have applied Equation [6] to the n(r) derived from measurements that were collected at 
fixed time intervals after the solution was made hypotonic. In Figure [6| these are compared with the 
normalised MSD of the same cell in phosphate-buffered saline (PBS) solution, prior to hypotonic exposure. 
Note that in order to more easily discern visually the changes in 11 (r), caused by the hypotonicity, the 
normalised MSD for the cell in PBS is plotted in each of the panels. 

The elastic {G'^{uj)) and viscous {G'^{uj)) components of the complex modulus, G'*{Ld), derived from 
the measurements shown in Figure |6] are shown in Figure [t] (scaled by the plateau modulus G^(0)). It 
is clear that, in PBS (Figure [7|\) , the high frequency behaviour of the cell's elastic and viscous moduli 
both show a frequency dependence of oc a;'^/'*, which is characteristic of isotropic in vitro reconstituted 
actin filament solutions |35l 136] . 

As the solution is made hypo-osmotic, the cell cytoskeleton starts to reorganise and the frequency 
behaviour of the moduli drastically change as (sequentially) shown in Figure ??. In particular, after 
~10min in the hypo-osmotic solution both moduli tend to assume a high frequency behaviour oc ui^^^ 
(see Figures [7]C and[7|l)). This change in the degree of elasticity may be explained by two alternative 
(and possibly concomitant) processes: (i) an increase in cytoskeletal tension as response to a stretching 
force caused by the cell swelling and (ii) the consequent formation of ^^dangling branches''^ [291 1301 131j 
(Figure |5| generated during the reorganisation process of the cell cytoskeleton (including those possibly 
obtained from the rupture of the cytoskeletal protein filaments because of stress overload). However, 
whilst the increase in cytoskeletal tension would explain both the initial increase of the cell stiffness 
and the change in the frequency scaling laws of the moduli (i.e. from cx Lo'^f'^ to cx o;^/^), it would fail 
to explain the softening process occurring during the cell volume re-equilibration, which ends with a 
breakdown, at high frequency, of the cell's elastic modulus (Figures [7]GI). This latter phenomenon could 
be simply justified by an excess of ^^dangling branches^\ Indeed, from a rheological point of view, these 
would only contribute to the capability of the cell to dissipate energy but not to the cell elasticity, since 
they cannot bear stress and hence do not contribute to the cell's rigidity. 

Finally, it is important to highlight that the results that have been obtained here using a video particle 
tracking method are in good agreement with those presented in the review written by Papakonstanti and 
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Stournaras [3], where a set of assessments of actin cytoskeleton dynamics and actin architecture in ceU 
volume regulation are summarised. However, none of the techniques reviewed in that work [3] are able 
to provide quantitative information on the cell viscoelasticity, as is the case here. 
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Figure 6. The normalised mean square displacement vs. lag time of a 5 iim diameter silica 
bead chemically bound to a Jurkat cell in iso osmotic (PBS) solution (square symbols) 
and in hypo-osmotic solution (circle symbols) after the addition of 10% v/v distilled water 
to the PBS buffer and measured at time intervals {At) of (A) 2 min, (B) 7 min, (C) 12 
min, (D) 18 min, (E) 23 min, (F) 28 min, respectively. At short time intervals, the n(r) has the 
potential to reveal the dynamics occurring at molecular level (as shown in Figure ??); whereas, at large 
lag-times it provides information on the stiffness of the whole cell (as shown in Figure |4]) . 
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Figure 7. The real (storage, G'^{uj)) and the imaginary (loss, G"(w)) parts of the complex 
modulus (G*(a;)) scaled by G'q vs. frequency of a Jurkat cell in iso-osmotic (PBS) solution 
(A) and in hypo-osmotic solution (B-G) after the addition of 10% v/v distilled water to 
the PBS buffer and measured at time intervals (A<) of (B) 2 min, (C) 7 min, (D) 12 min, 
(E) 18 min, (F) 23 min, (G) 28 min, respectively. The moduli have been evaluated by using 
Equation [6] on the normalised mean-square displacement data shown in Figure |6] The lines are guides 
for the gradients. Note that, the scatter in the data at high frequencies and the ripples in the low 
frequency portions of the curves in Figure [t] are due to the analytical method that we have adopted |19| 
for performing the Fourier transforms of the normalised MSDs. As explained in detail in Ref. [37], this 
method works directly on the experimental data points (i.e. |rfe, life}, where k — 1..N) and therefore 
preserves genuine experimental noise. 
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Conclusions 

In summary, we have presented a straightforward and non-invasive experimental procedure, coupled 
with a new analytical method to interpret the data, which leads to quantitative determination of the in 
vivo viscoelastic properties of cells in the frequency domain. The method has the potential to monitor 
the internal dynamics and re-organisation of the actin cytoskeleton up to frequencies on the order of 
kHz, representing a valuable addition to studies that address cellular physiology and pharmacological 
response. Indeed, in this work we report, for the first time, the high frequencies (up to ~600Hz) changes 
of the Jurkat cells' viscoelastic spectrum from oc w^/^ to oc w^/^, as response to a change in osmolarity 
of the solution. The rhcological interpretation of the results gives a direct insight of the coll cytoskeleton 
structure and its re organisation. In the future, it is envisaged that these interpretations could be 
coupled with advanced molecular biology techniques to resolve the detailed interactions underlying these 
rheological changes and that faster dynamics could be studied by means of a quadrant photo-diode based 
tracking system 
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